デバ イス の 記事 


工 
⑧ 


今回 は , 乗算 を ベー ス に する 潤 近 法 を 用 いた 除算 回 路 の 実現 方 
法 を 解説 する . 具体 的 に は , 除算 に 特 化し た 「 反 復 乗 算法 」 と 
汎用 的 な 「Newton-Raphson 法 」 の 二 つ の アル ゴリ ズム を 用 
いた 回 路 設 計 を 紹介 する . な お , 本 稿 で 設計 し た 回 路 の ソー 
ス ・ コ ー ド は 本 誌 ホ ー ム ペー ジ (http://www.cqpub.co.jp/ 
dwm/) か ら ダ ウン ロー ド で きる . (編集 部 ) 


これ まで 3 回 に わた っ て , 回 復 法 や 非 回 復 法 を 用 いた 除 
算 回 路 設計 お よび 高速 化 に つい て 話 を し て きま し た . これ 
ら は , いずれ ぞ 被 除 数 の 中 に 除 数 が いく つ 取 れる か 」 と い 
う 視点 に 立っ た も の で し た . 
今回 は ,「 被 除 数 と 除 数 の 比 を 求め る 」 と いう 視点 で 除算 
実行 する 方 法 に つい て 考え て み ま し ょ う . 代表 例 と し て , 
激 近 法 を 取り 上 げ ま す . 


[1 新法 と は 


た し か 高校 の 数 学 だ っ た と 思い ます が , 尊 化 式 と いう も 
の を 習っ た と 思い ます . 初期 値 ,。 および, 番目 と 証 1 番目 
の 変数 の 関係 を 定義 する こと で , すべ て の 変数 を 定義 で き 
る と いう も の で す . ヵ - ゅ oo と し た と き に 収束 する か , な ど 
も 議論 され て いた よう に 記憶 し て いま す . 

潤 近 法 は この 考え か た を 応用 し た も の で す . 初期 値 か ら 
は じ め て 番目 か ら j+ 1 番目 の 変数 を 求め る こと を 繰り 返 
し , 特定 の 値 に 収束 させ る こと で , 関数 の 近似 を 行う と い 
う も の で ず も ちろ ん , 無限 に サイ クル を 繰り 返す の は 現 
実 的 で は な いた め , 有限 回 の サイ クル 後に 必要 な 精度 を 得 
る も の と する ). 洒 近 法 の 中 に は , 例え ば Newton-Raphson 
注 詳細 は 後述 ) に 代表 され る 汎用 性 の 高い アル ゴリ ズム も 
あり ます が , 適用 する 関数 に よっ て く ふう が 活 ら さ れ た も 
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新 近 法 を 用 いた 陰 距 回 中 


関連 デー タ 


の も あり ます . まず は , 除算 に 特 化 し た アル ゴリ ズム を 紹 
介し まし ょ よう 。 


人 @ 回 復 法 も 非 回 復 法 も 実は 滞 近 法 の 一 種 
と ころ で , 回 復 法 を 用 いた 除算 回 路 で は 次 の 式 を 繰り 返 
し 演算 し て いた こと を 思い 出し て くだ さい . 


pg, 三 1 
の の ド py 生 の の 場合 ) で で orei…i… {( 1-1) 
am ニー 
pg, =0 

( く g の 場合 ) で Oren 《 1.2) 
7 な 1 モア 


ここ で , : 部 分 余 , pg,: 部 分 商 ,7: 除 数 
実は , 条件 分 岐 が あめ る と は いえ , 回 復 法 や 非 回 復 法 な ど 
も 洒 近 法 の 一 種 な の で す . けた 単位 で の 潤 近 を 反復 する こ 
と で , 商 を 必要 と する 精度 に 追い 込ん で いき ます . し か し , 
手段 に 筆算 の イメ ー ジ が 強い た めか , ある い は 結果 が 商 と 
余 で 誤差 な し で 与え られ る た めか , 潤 近 法 と し て 語ら れる 
こと は な いよ う で す . 


人 @ 乗算 ベー ス の 洛 近 法 
回 復 法 や 非 回 復 法 は , けた 単位 で 加算 減算 ) を 行う 方 法 
で し た . つま り , サイ クル 当たり 1 ビッ ト の 商 を 求め よう 
と すれ ば , た ビッ ト の 結果 を 求め る の に た サイ クル 必要 と な 
る の で す . 

これ に 対し , 今回 は 乗算 を ベー ス に する こと で 指数 的 に 
精度 を 上 げ ら れる 方 法 を 考え て み ま し ょ う . ここ で は ,「 反 
復 乗算 流 ある い は , Goldschmidt の アル ゴリ ズム )」 と 言 
われ る 方 法 を 用 いる こと で , 反復 する サイ クル を logzx に 
短縮 で きる こと を 説明 し ます . 

商 を 2, 被 除 数 を z, 除 数 を 7 と する と , 7 三 z/ の 9 と な り 


ト に に て @ て 
ニーー ーー ーー 一 デ ドー 
て て と に に 介 て に @ 介 
ーー 一 | 一 王 ーー トー コ に ーーー 
て て て ト - に 


ます . つま り , 求め る も の は z と 2 の 比 と な り , 式 2 の 表 
現 を 用 いま す . 


まず , この 式 が 何 を 表し て いる か を 説明 し ます . これ は , 
「: サ イク ル 目 で xx と いう 変数 を 分 母 と 分 子 そ れ ぞ れ に 乗算 
する 」 と いう 操作 を 反復 する も の と 考え て くだ さい . 分 母 
と 分子 に 同一 の も の を 乗算 し て いる の で , 2 と し て は 変化 
し な いこ と に 着目 し まし ょ う . ここ で , j サ イク ル 目 に お 
ける 分 母 を み , 分 子 を z, と する と , 式 3) が 成立 し ます . 


の の ・ 脳 
es 半 | ( 3) 


これ が 洒 化 式 と な り , 初期 値 は 2 と z で す . を 乗算 す 
る 操作 を 反復 する こと で みつ 1 と いう 滞 近 操作 が で きれ ば ぱ ば, 
つの と いう 結果 を 得る こと に な り ま す . つま り , 分 母 を 1 
に 追い 込ん で いく こと で , 自動 的 に 分 子 が 2 に 追い 込ま れ 
る と いう わけ で す . 

で は , を どの よう に 定義 むせ する か を 説明 し まし ょ う . ま 
ず , みっ 1 を 実現 する こと は 1- みっ 0 を 実現 する こと と 等 
価 と 考え る こと が で きま す . 式 2) の 分 子 の 初期 値 は , 


1/2 ミ 2 ぐ 1 CPPPPPPPFFFF ど FFP です で すす どど すす す と と すす すす ミド すす すす で ( 4) 


に 正規 化し て お きま す . する と , 


0< く 1- た 1/2 PPPPPPPFPFPFPF で PP で PP ど た と と どど と すす ミ すす ミミ すす すす ( 5) 


と な り ま す . この た め , 1- 97 の べき 乗 で あれ ば 0 に 近づけ 
る こと が 可能 で ある と 判断 で きま す . ここ で は 乗算 を 用い 
て 洛 近 操作 を 行い ます が , 次 数 を 高く と る と 乗算 回 数 が 多 
く な り ま す . そこ で , 次 数 は 2 次 乗算 1 回 ) と し て , 以下 
の よう に 治 化 式 を 定義 し まし ょ う . 


7/EE 拓 ) 7/5 REEEERRREECOOEEEOEEREEEEEEEREEEEE ( 6) 
式 6) を 変形 する と , 
1- の] ネ | 6MRSGNASS2D8GReRg 《 の ) 


と な り , サイ クル が 進む ほど 1- みつ 0 を 実現 で きる こと 
が 確認 され まし た . また , 式 6) は 次 式 の よう に 変形 で き 
ます . 


人 CC イィ ン ン ーー 4 
0 っ = な ご ーー 
の (7) デー に /) 町 
し ー (62P 8) 


葉 8) に 3) を 代入 する と , 欧 9) が 得 ら れ ま す . 


これ を 順次 み , に 乗算 し て いく の AG みつ 1 る の 
が 実現 で きま す . 


人 @ 諾 差 が ある こと を 考慮 し て 変数 の ビッ ト 幅 を 決め る 
演算 回 路 ハー ド ウェ ア ) に お ける 変数 の ビッ ト 幅 を と 
する と , み は 1- 2 * 以 上 に は 1 に 近づく こと が で きま せ 
ん . し た が っ て , 2 三 た を 満た す ,7 に 至っ た 時 点 で 演算 の 反 
復 動作 は 終了 と な り ま す . この こと か ら , 反復 回 数 : は , 次 
の 式 で 与え られ ます . 


7 三 og ん を 満たす 最 小 の 整数 。…‥……………|… ツ 。 10) 


回 復 法 や 非 回 復 法 な どの 場合 , 除 数 , 被 除 数 , 商 , 余 の 
関係 に は 誤差 が あり ませ ん で し た. 一 方 , 反復 乗算 法 で は 
「 得 ら れ た 4 商 ) は 基本 的 に 誤差 を 含ん で いる 」 と いう こと 
に 注意 し な く て は な り ま せん . 演算 回 路 に 誤差 が な か っ た 
と し て も , 演算 方 法 そ の も の に 「 み は 1- 2 * 以 上 に は 1 に 
近づく こと が で き な い 」 と いう 大 前 提 が ある の で , 分 子 の 2 
も その 分 小さ な 値 を と る こと に な る か ら で す . 

また , 乗算 を 行っ た 結果 , 変数 の ビッ ト 幅 は 多 と な り ま 
す . これ を 次 の サイ クル に その まま 持ち 込む と , 変数 の 幅 
が サイ クル を 重ね る た びに 大 きく な っ て し まい , 現実 的 で 
は あり ませ ん . 乗算 結果 は ビッ ト に 丸め て , 次 の サイ ク 
ル に 渡す の が 一 般 的 で し ょ う . この 丸め に よっ て も 誤差 が 
生じ る こと に な り ま す . 丸め に よる 誤差 を = と する と , 誤 
差 を 含ん だ は 式 11) の よう に 表現 で きま す . 


gg 三 の 十 りつ 間 還 UCCCLLCLKK エ ミミ ミエ ナミ ミエ ミミ エト ミ ミエ ナミ ミミ ミミ エミ ミミ ミミ ミミ ミミ ミミ で ( 11) 


蘭 11) を 用 いて 式 3) を 書き 直す と , 葉 12) が 得 ら れ 
ます . 


みみ EX 2ー カー と ) 

三 w 2- の 十 28 1- の 填 2 
Zr1 王 人 4 2 み ー と ) 

= テス 2- 忌 十 gz, 


PC ( 12) 


丸め を 単純 な 切り 捨て で 行っ た と する と , = ぐ 1ul unit 
mm least position : 最 下位 ビッ ト ) と な り ま す . =2 は 十分 小 
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さく 無視 で きる も の と 考え る と , 誤差 と し て 支配 的 な の は 
葉 12) の 灰色 の 部 分 と な り ま す . この 誤差 は サイ クル を 引 
ね る た びに 累積 され る た め , 求め る 商 の けた に よっ て 反復 


Im 
IIU 


[2 反復 乗算 流 の 回 路 化 


回 数 を 決定 し た 後 , 累積 誤差 の 最大 値 を 考慮 し て 演算 回 路 1 に , 反復 乗算 法 に よる 除算 器 の ブロ ッ ク 図 を , リ スト 1 
に お ける 変数 の ビッ ト 幅 , を 決め る 必要 が あり ます . に その Verilog HDL の RTL 記述 を 示し ます . 
d[232:0] z [23:0] 図 4 
邊 
dr[31:0] 人 
尽 
王 
31 24 16 8 0 31 24 16 8 0 
dr[31:0] ご zr[31:0] 
当 1 z [31:0] 
2- di 図 
半 
1 1 
い 


*_mul[63:0] 図 
( 乗算 結果 ) 


32 ビ ッ ト ( zr, dr に フィ ー ド バッ ク す る 範囲 ) 較 


| = 


zr[31:0] 


人 箇 生 作 | 


q[31:0] 


M 


峰 


1 反復 乗算 法 に よる 除算 器 ブロ ッ ク 
に は , 前 処理 . ルー ズ 反復 ), 後 処 理 と ある が , 実際 に 演算 処理 を 行っ て いる の は ルー プ 部 の み . 


較 


リス ト 1 図 1 の 回 路 構 成 の RTL 記 述 div_convergence_a.v) 


modu1e diy oonvergenoe a ( if ( re ) ュ <= 3'hO 
d // quotienE [23 : 00] : Eraction elge FE ( gtart ) 1i <= 3!h5 
Dusy // busy else 1E ( en sed ) 1 <= 1 - 1'b1 
end 


7 


diyidend [23 : 00] : Eraction 
diyisor [23:00] :Fraction // Eunction (演算 処理 部 
Ear 上 (Qivider 上 rtigger ) wire [31:00] d cmp1 cdqr + 1"b1  : // gr の 2 の 補 数 
clock wire [63:00] d mult dr * d mp] 』 
エ reset ena // 分 母 に qr の >2 の 補 数 を 乗算 
wire [63:00] 2 mult Zr * d omDp] : 

// 分 子 に qr の 2 の 補 数 を 乗算 
: 00] H // nexxt reg contents (次 の サイ クル の た め の レ ジス タ 内 容 ) 
: 00] : wire [31:00] dw d mu1t[62:31] : // 分 母 の 丸め 
wire [31:00] zw z_mu1t[62:31] : // 分 子 の 丸め 


a1way8 @( posedge cl1k ) begin 
1E ( start ) zr <= { 1'b0,z,7!h00 】 : 
e1ge 1F ( en sed ) zr <= zW 7 


end 


a1way8 @( posedge clk ) begin 
[02:00] : H ( start ) dr <= { 1'b0,d,7'h00 】 : 
[31: 00] : // [31] :int [30:00] :Erac else 1E ( en seq ) dr <= dw : 
[31:00] : // [31] :int [30:00] :Erac end 


! 


// sequence (シー ケン サ ) agg1gn d gr [31:08] + zr[07] : // 分 子 を 丸め て 出力 


wtre en seqd= (1>0) : a8819n DuSY en Sed 


a1way8 @( posedge olk ) begin endmodu1e 
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と て て ト て 
ーーー ーー ーーー 一 | 王 デ トド ーー ーー 
て 1 し 1 に 介 に に 6 介 
ーー 一 | 一 ーー ーー トー し ーーー 
と に て @ し 


ここ で は 式 4) を 満た す た め の 入力 の 正規 化 は 済ん で い 
る も の と し ます . 浮動 小数 点 演算 の 様相 を 呈し て いる ので, 
z に つい て も 同じ よう に 正規 化 を 施す こと と し , 単 精度 の 
演算 を 意識 し て 入力 の ビッ ト 幅 を 24 ビ ッ ト と し まし た . 式 
( 10) よ り 反 復 回 数 は 5 と な り ま す . 反復 回 数 を 5 と する と 
32 ビ ッ ト まで 演算 で きる の で , 内 部 演算 の ビッ ト 幅 は 32 
ビッ ト と し まし た が , 実際 に は 28 ビッ ト 程度 で 十分 な は ず 
で まま 

ブロ ッ ク 図 は , 前 処理 , ルー ズ 反復 ), 後 処理 に 分 け て 
描い て あり ます が , 前 処理 は 24 ビ ッ ト の 入力 の 32 ビッ ト ・ 
レジ スタ へ の 割り 当て を , 後 処理 は 32 ビ ッ ト ・ レ ジス タ か 
ら 24 ビ ッ ト の 商 を 取り 出す 位置 を 示し て いる だ け で す . す 
べ て 反復 動作 の み で 除算 が 完了 し ます . 


@ ある 程度 の 精度 で あれ ば FPGA の ハー ド ・ マ クロ を 利用 
乗算 器 は ター ゲッ ト ・ デ バイ ス に よっ て 実現 方 法 が 異な 
る た め , リス ト 1 の RTL 記述 で は 単純 に 乗算 を 示す 演算 子 
の み で 表現 し て いま す . また , 乗算 は 1 クロ ッ ク で 処理 で 
きる こと を 前 提 と し て いま す . この まま で も ター ゲッ ト ・ 
デバ イス に よっ て は 実用 可能 な は ず で す が , あく まで も 演 
算 方 式 を 理解 し て も ら う た め の 参 考 記 述 と お 考え くだ さい . 
シミ ュ レ ーション で 分 母 と 分 子 の 尊 近 動 作 を 確認 し て も ら 
えれ ば , と 思い ます . 

乗算 反復 法 を 実際 に 用 いる 場合 , 動作 速度 を 決定 する の 
は 2 の 2 の 補 数 演算 と 乗算 の 部 分 に な り ま す . ここ が 設計 
者 の 腕 の 見 せ ど こ ろ と な る で し ょ う . シー ケン ス の くさ ふう 
で 乗算 器 を 1 個 で 済ま す の も お も し ろ い と 思い ま す . 

図 2 に , シミ ュ レ ーション 結果 の 例 を 示し ます . 7= テ 05, 
z 三 0.75 の 場合 , 商 は 1.5 に な り ま す . 反復 回 数 が 大 きく な 
る ほど 分 長 の ⑦) は 1 に , 分 来 >) は 1.5 に 近づい て いる の が 
わか り ま す . 

乗算 反復 法 は 乗算 器 を 用 いる と いう こと で , 大 げ さ な 方 
式 と 思わ れ が ち で す が , 近ごろ の FPGA の よう に 乗算 器 が 
ハー ド ・ マ クロ と し て 用 意 さ れ て いる の で あれ ば , 比較 的 
採用 し や すい の で は な いか と 筆者 は 考え ます . も っ と も , 
高 精 度 な も の が 必要 と な る と , 乗算 紅 も 大 き な も の が 必要 
に な る の で , 適用 範囲 に は 限界 が ある の か も し れ ま せん . 


[3 Newton-Raphson 法 に よる 逆数 演算 


本 稿 の 最初 に , 潤 近 法 の アル ゴリ ズム の 例 と し て Newton- 


。 凡 ら 健 回 褒 諾 針 。 


に ら 一 ー 三 シシ ジ m 三 
(0( 品 = ご 0⑰ 記 = 
Raphson 法 の 名 を 挙げ まし た . この 方 法 は 著名 な の で す が , 


除算 に 適用 し よう と する と めん どう で す . 除算 は 引き 数 が 除 
数 と 被 除 数 の 二 つ で す が , Newton-Raphson 法 は 微分 を 使用 
する た め , 2 変数 で は ぐ あい が 悪い の で ず 不可 能 で は な い 
が , 除算 回 路 の 実現 と いう 目的 か ら は 適当 で な いと いえ る ). 

し か し , 除算 を 被 除 数 と 除 数 の 逆数 の 乗算 」 と と ら え れ 
ば , 除 数 の 逆数 を 求め る の に Newton-Raphson 法 は 有効 で 
す .「 逆数 を 求め て , さら に 乗算 する な ん て 二 度 手間 だ 」 と 
思わ れる か も し れ ま せん が , 例え ば アフ ァ イ ン 変換 に よる 
座標 変換 な ど で は 共通 の 分 母 で 複数 の 分 子 を 除算 する と い 
っ た 操作 が 行わ れ ま す . この よう な 場合 , めん どう な 除算 
は 逆数 と し て 1 回 だ け 求 め , 後 は それ ぞ れ の 分 子 に 乗算 す 
る と いう 方 法 が 有効 で す . 


@ の 接線 と x 軸 の 交点 で 滞 化 式 を 定義 する 

それ で は , 逆数 演算 に Newton-Raphson 法 を 適用 し て み 
まし ょ う . まず は アル ゴリ ズム か ら で す . 図 3 に その 原理 
を 示し ます . ここ で , 式 13) に 示す 関数 を 定義 し ます . 


た ェ ) 三 1/*ー 人 CCCCCCCPPPPPPPPPPPHTPEPPPPEPPPPPPPPPPHPPEP ( 13) 


式 13) に お いて , fx>) テ 0 と お く と ェ x テ 1/2 を 根 に 持ち 
ます . 言い 換え れ ば , x+= 1/7 で x 軸 と の 交点 を 持ち ます . 
この よう す が 図 3 に 示さ れ て いま す . 

さて , + テ ニャ ェ に お ける テル た >) へ の 接線 を 定義 する と , 
図 3 に お ける 直線 ヵ に な り ま す . 直線 7 ヵ と x 軸 の 交点 を 
xn1 と する と , テテ ェ ュ に お ける 新た な 接線 な が 定義 され , 
これ と x 軸 の 新た な 交点 x+2 が 求め られ ます . 図 か ら も わ 
か る よう に , この 操作 を 繰り 返す こと で が 求め られ る ベ 


6 
14 
1.2 
1.0 ー ャ ーーーー 
0.8 了 還 記 計 
0.6 6 
』 ブ 
0.4 
0 1 2 8 4 5 


反復 回 数 較 
図 2 乗算 反復 法 に よる 分 母 と 分 子 の 滞 近 


g 三 0.5, > 0.75 の 場合 , 商 は 1.5 に な る . 反復 回 数 が 大 きく な る ほど 分 母 
( 9) は 1 に , 分 求 z) は 1.5 に 近づい て いく の が わか る . 
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て 


ーーー ニー 
に 


と 1 に 生 得 G て に に に 
1 1 に ( ( 信 に 1 人 介 
C < C C C て に て 6 @ G C @ C @ て で @ @ @ 
い W し 1 に し 【 6 1 W 
て と し 旭 と 6 し 休 て @ に 
ーg 1/ ガ ー ェ ) 2 
ピア き 20 じ POPHHEPEPPPPPREEEEEEEPPPPEEERPRPEEEE ( 18) 


2 7+1 +2 


図 3 Newton-Raphson 法 に よる 逆数 演算 の 原理 
関数 f x) ニ 1x- g は , 1 が x 軸 と の 交点 と な る . 


き ェ = テ 1/7 に 近づい て いき ます . つま り , と 1 の 関係 を 
定義 する 潤 化 式 と , 初期 値 xy が 定義 で きれ ば , 滞 近 法 と し 
て 逆数 演算 が 可能 に な が り ま す . この よう に に お ける 接線 
と x 軸 の 交点 で 潤 化 式 を 定義 する の が Newton-Raphson 法 
で す . 

で は , 潤 化 式 を 求め る こと に し まし ょ う . まず , x テ ェ 
に お ける ッ y テ ん / >) へ の 接線 は , 式 14) で 定義 され ます . 
放 ) は * デ ニャ に お ける 接線 の 傾き を 導 関 数 を 用 いて 表現 し 
た も の で す . 


ES あう た 4 の うま だ まう 間 ES 《 14) 


この 直線 と x 軸 の 交点 x+1 は 下 14) に お いて = テ 0 と し 
た と き の 解 な の で , 式 15) の よう に 表せ ます . 


re aa ( 15) 
ここ で, 
誠 ェ ) ーー 1/x 3 858PS3S0SISiS88iEIS お eS133829398IG13208197818feRSOEI8288 が な ( 16) 


な の で , 式 15) に 式 13) と 葉 16) を 代入 する と , 次 の 手 
化 式 が 得 ら れ ま す . 


621 本 PERPP 生 は 相 P ( 17) 


後 は , 初期 値 が 定義 で きれ ば 逆数 演算 が 成立 し ます . 
ヵ ス テッ プ 目 の 誤差 を e, 三 1/2- xx と する と , 


e+ ュー 1/ ガ ー 4 ュ 
三 1/29- XX 2- ェ ・ の 
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ここ で , 反復 乗算 法 と 同じ く , 式 4) の よう に 正規 化し 
て お く MG gr1 く ( g,) 2 と な り ま す . 式 4) に 従っ て , 


1 ぐ xo ミ ウ OCCCCCKKTKK エ エメ エエ ミエ エミ エエ ミエ エミ エエ ナ ミエ ミミ ミミ ミミ ミド すま すす ます ( 19) 


と な る よう に 選択 され て いれ ば , 洒 近 が 保証 され ます . つ 
まり , 初期 値 は 7 の 値 に 関係 な く , 例え ば 1.5 の 固定 値 と し 
て も か まわ な いわ け で す . 

し か し , 式 18) か ら 明らか な こと で す が , 7 ステ ッ プ 目 
に お ける 商 の 精度 を , ビ ッ ト と する と , :+ 1 ステ ッ プ 目 で 
は 2z, ビ ッ ト 以上 の 精度 と な る の で , x+o の 選択 は その 後 の 
反復 回 数 を 左右 する こと に な り , 重要 で す . 

例え ば , 4 ビッ ト 精度 の xo で 演算 を 開始 すれ ば , 3 回 の 
反復 で 32 ビ ッ ト 精度 まで 追い 込む こと が 可能 と な り ま す . 
以上 か ら , 所 望 の 商 の ビッ ト 精度 を ん , 初期 値 の ビッ ト 精 
度 を zo。 と する と , 必要 な 反復 回 数 : は 次 式 で 求め られ ます. 


7 且 log6 4/go) を 満た す 最 小 の 整数 …………… ( 20) 


⑱ テー ブル を 使っ つて 適切 な 初期 値 を 与え る 

4 に Newton-Raphson 法 に よる 逆数 器 の ブロ ッ ク 図 を , 
リス ト 2 に その Verilog HDL の RTL 記述 を 示し ます . 反 
復 乗算 法 の 場合 と 同じ く , 式 を 満た す た め の 入力 の 正 
規 化 は 済ん で いる も の と し , 入力 の 条件 も 24 ビ ッ ト 幅 と し 
まし だ ブロ ッ ク 図 お よび RTL 記述 に お ける 注意 点 も 反復 
乗算 法 の 場合 と 同じ ). 

すでに 説明 し た と お り , Newton-Raphson 法 で は 初期 値 
を より 良い 精度 で 与え る こと で 反復 回 数 を 少な く で きま す . 
ここ で は , 1.5 な どの 固定 値 で は な く , テー ブル を 用 いて 初 
期 値 を 適切 に 与え る こと に し まし ょ う . 

し か し , テー ブル が あま り に か さば り すぎ る の も 問題 で 
す . ここ で は 4 ビッ ト 精度 で 初期 値 x ぃ を 与 る る こと に し ま 
す . リス ト 2 の Eunction 文 に ある case 文 が こ の テー ブ 
ル を 表し て いま す . 1/2 ミ 2 ぐ 1 の 範囲 を 24 ビ ッ ト 表現 で 
入力 する の で , MSR most significant bit) は つね に 1 と 
な り ま す . し た が っ て , 』 22: 19] で テー ブル を 索引 すれ 
ば よい こと に な り ま す . 本 来 , 出力 の 範囲 は 欧 19) の と お 
り で す が , 1 ぐ xop く 2 と な る よう に テー ブル の 内 容 を 選択 
すれ ば , や は り MSB は つね に 1 と な り ま す . この た め , 


| に に / ュ ィ っ 選ん 
ce C 6 C C _@E 、 角 9 企 三 回 有 E 半 る ia 
C G C @ e@ @ C C C C ー 4 
ーー 一 | ーー ーーー トー コ に ーーー: プ 
し ーー (-、 一 ご 一 
@ C C @ C (/ し ーー (7 ミー 
の (02=C_ OO 弓 = 
d[23:0] 1 
dr[22:19] 初期 値 テ ー ブ ル 図 較 
( 4 ビッ ト x16) 図 に 
軸 
小数 点 位置 拓 
尽 
dr[23:0] 
23 16 8 0 31 24 16 8 0 1 
1 1 
31 24 16 8 0 
zr[31:0 
32 ビ ッ ト ( zd_cmp で 2 の 補 数 を と る 範囲 
較 
6 zd[55:0] 生 
zd[55:0] 図 
( 乗算 結果 ) 較 2- zd 貞 
zd_cmp [31:0] 
zw [63:0O] 図 
( MM ー 32 ビ ッ ト ( zr に フィ ー ド バッ ク す る 範囲 ) 凶 J 
員 56 48 40 32 
4 
年 
zr[31:0] 回 
24 島 
q[23:0] 燥 


図 4 Newton-Raphson 法 に よる 逆数 器 ブ ロッ ク 
反復 回 数 を 少な くす る た め 


リス ト 2 5 の 回 路 構 成 の RTL 記述 div_convergence_b.v) 
modu1e diy convergenoe_D ( 
d , // quotient [23:22] 
bugy // busy 


:1nE [21:00] : Fraction 


// divigsor [23:00] :Fraction 
/ / star 上 (divider 上 て 1qqder) 
// c1ock 

: // rese ena 


[23:00] d 
Star : 


Cc1k 


[02 : 
[31: 
[23: 


00] 1 : 
00] zr : // 
00] gr : // 


: エエ aC 
: エエ aC 


[31:30] :int [29: 


[23 : 


00] 
00] 


(シー ケン サ ) 
(>0): 


// sequenoe 
wire en SeQq = 


a1way8 @( posedge Clk ) 
ュ E ( re ) 1 <= 3'hO 
e1ge 1 ( gar ) < 3Th3 
e1ge FE ( en seqd ) 1 <= 1 - 1!D1 : 


begin 


end 


// Eunction (演算 処理 部 ) 
// 初期 値 テ ー ブ ル 
Eunoction [03:00] finix0 : 
nput [22:19] d uppr : 
Ca8@e ( d_uppr ) 
4「!h0 : Finix0 
4!h1 : finix0 = 


4'hE 
4'hD : 


め , ここ で は 初期 値 と し て 固定 値 で は な く , テー 


ブル を 


用 いて 適 


|= 


な 値 を 与え る こと に し た . 


4!h2 : finmix0 = 
4「h3 : finix0 
4「!h4 : finix0 
4!h5 : finix0 
4!h6 : fimnix0 
4「!h7 : finix0 
4!h8 : finix0 
4!h9 : finix0 
4!ha : finix0 
4!hb : finix0 
4!ho : finix0 
4!hd : finix0 
4!he : finix0 
4「!hfF : FinixO = 
endoag8e 
endFunot1on 
wire [29:26] z0 = Finix0 (d[22:19] ): // 初期 値 の 決定 
// 1oop 
wire [55:00] zd = アア * dP 
// xi*d [55:54] 
zd Cmp = “zd[54:23] 
// 2-xix*d [31] : in モ 
wire [63:00] zw = Zr * zd Cmp 月 
// xi*(2-xt*d) [63:61] : inE 


: エエ ao 


7 
Re eoN 
+ 1!b1 : 


[53 :00] 
wire [31:00] 
[30 : 00] : Erac 


[60 : 00] : Erac 
a1way8 @( posedge clk ) begin 
ュ E ( start ) gzr <= { 2'h1,z0,26'h0000000 }: 
e1se 1fF ( en seq ) 2 <= 2W[62:31] : 
end 


a1way8 @( posedge c1k ) begin 
1E ( gEarE .) qr <=dz 
end 


= [31:08] 
en sed : 


a881gn d 
a881qn Du8Y 


+ Zr エ [07] : 


endmodu1e 
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と て 1G 生 1 の て て 6 に 
し に 《 ( 細 ( て = し 1 
て て に 1 に 介 に 6 6 介 て トレ て 介 で で て に 6 
ーー | 一 【 1 @ 1 ( 1 も 6 
C 【@ て @ ( に で C @ @ 
2 
凡人 
較 
回 
宰 
屋 1.5 
SM ーー 初期 値 4 ビ ッ ト 精度 図 
ーーー 初期 値 1.5 固 定 較 H 
1.25 1.4 1 1 
0 1 2 3 4 
6 | 上 本 国人 
0.5 0.625 0.75 0.875 1 図 6 初期 値 の 違い に よる 滞 近 の よう す 
す 23:19] 図 24 ビ ッ ト 精度 の 商 を 求め て いる の で 図 か ら は わか り に く い が , 初期 値 を 固 


5 逆数 演算 初期 値 テ ー ブ ル 

4 ビッ ト 幅 で 構成 し た 場合 の MSB を 含む テー ブル の 内 容 を 示す . な お , こ 
の テー ブル の 内 容 だ と , 4 ビッ ト 精度 に は 若干 不足 する 場合 が ある の だ が , 
実 使用 時 に 問題 が な か い の で この まま と し て ある . 


この 部 分 は テー ブル に 含ま ず ザ レ ジス タ の 初期 化 で 与え る こ 
と と し , 4 ビッ ト 幅 で テー ブル の 内 容 を 構成 し ます . MSB 
を 含ん だ テー ブル の 内 容 は 図 5 の よう に な り ま す . 

以上 の 条件 か ら , 式 20) に より 反復 回 数 は 3 と な り ま す . 
また , 図 6 に 初期 値 を 変化 させ た 場合 の 滞 近 の よう す を 示 
しじま す 。 


今回 は , 潤 近 法 に 焦点 を あて て み ま し た . 設計 者 と し て 
駆け 出し の ころ , 筆者 に と っ て も この よう な 乗算 の 手法 は 
と っ つき に くく , で きる だ け 人 避 ける よう に し て いた 記憶 が 
あり ます . し か し , 数 式 に 惑わ され ず , ツボ さえ 押さ えて 
お け ば めん どう な も の で は な いこ と が 理解 で きた こと と 思 


定 の 1.5 と する と 4 ビッ ト 精度 の テー ブル を 用 いた 場合 より も 反復 回 数 が 一 
つ 多 く な り , 4 回 と な る . 


いま す . 精度 を 確保 する 反復 回 数 も , 意外 と 単純 な 計算 で 
求め られ ます . 構造 的 な 面 で 特殊 な の で , オー ル マ イ ティ 
と は 言い が た い の で す が , 条件 が 合え ば 使っ て みる の も お 
も し ろ い で し ょ う . 


すず き ・ し ょ うじ 
スタ ンダ ー ド ・ ロ ジッ パル 株 ) 


筆者 プロ フィ ー ル ツ 

鈴木 昌治 . 10 月 初旬 に ,「 グラ ン フ ォ ン ド 終 魚 川 120km ほど の 山 
あり 谷 あ り コ ー ス を 走破 する 自転 車 遠足 ) に 参加 . ひど い 雨 の 中 で の 
走行 と な り ま し た が , 参加 賞 と し て カニ を も ら っ て ホク ホク . 
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右 の 書 に な る こと で し ょ う . 


ハー ド ウェア 記 述 言語 の 速 習 & 実 貴 


改訂 入門 Verilog HDL 記述 


小林 優 著 B5 変 型 判 256 ペー ジ 定価 3,.360 円 税込 ) ISBN4-7898-3398-4 
本 書 は , 回 路 図 ベー ス の 設計 か ら HDL に よる トッ プ ダ ウン 設計 に 移行 し た いと いう 方 , ある い は 初め て HDL に よる LSI 設 
計 に 携わる 方 な ど に 最適 如 Verilog HDL」 の 入門 書 で す . 実例 が 豊富 に 掲載 され て お り , FPGA や ASIC を 設計 する と き の 座 
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